Crossed Andreev Reflection in Quantum Wires with Strong Spin-Orbit Interaction 
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We theoretically study tunneling of Cooper pairs from an s-wave superconductor into two semi- 
conductor quantum wires with strong spin-orbit interaction under magnetic field, which approximate 
helical Luttinger liquids. The entanglement of electrons within a Cooper pair can be detected by the 
electric current cross correlations in the wires. By controlling the relative orientation of the wires, 
either lithographically or mechanically, on the substrate, the current correlations can be tuned, as 
dictated by the initial spin entanglement. This proposal of a spin-to-charge readout of quantum 
correlations is alternative to a recently proposed utilization of the quantum spin Hall insulator. 
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One of the key features and resources of quantum 
mechanics is entanglement, particularly in the particle 
spin sector, which has been an enticing subject since 
the Einstein-Podolsky- Rosen thought experiment [T] and, 
more recently, fueled by the modern proposals for spin- 
based quantum information processing and computa- 
tion [2]. In order to use an entangled pair of elec- 
trons for quantum information technology in a scalable 
semiconductor setting, it is essential to have a solid- 
state system that can separate the entangled electrons 
over appreciable distance. Detecting electron spin en- 
tanglement is possible via bunching or antibunching cor- 
relations in beam splitters [5] and transport through 
Coulomb-blockaded quantum dots forming a Josephson 
junction 0]. A conceptual headway came with a proposal 
to spatially separate spin-singlet Cooper pairs (CP's) in- 
jected from an s-wave superconductor via crossed An- 
dreev reflection (CAR) 5J in a normal- metal fork [5J. 
Later, more elaborate considerations for an s-wave su- 
perconductor in junction with quantum wires |S] and 
quantum dots [5] have been put forward. CAR is es- 
sential in all these proposals, and it has been experi- 
mentally manifested in the negative nonlocal differential 
resistance in the system of superconductor in junction 
with normal metal 10J. CP splitter experiments have 
been recently performed with quantum dots and car- 
bon nanotubes [123 • ^ s another form of CP splitter, we 
theoretically proposed the system with superconductor 
straddling a strip of two-dimensional quantum spin Hall 
insulator (QSHI) (T3], to inject a CP into its gapless edge 
states. Utilizing the helical Luttinger-liquid (LL) charac- 
ter of the QSHI edges (where each electron moves in the 
opposite direction to its time-reversed Kramers partner 
with opposite spin), the spin entanglement can be con- 
verted into nonlocal charge-current cross correlations. 

In this paper, we consider CP injection into quantum 
wires with strong spin-orbit interaction (SOI), such as 
self-doped (and possibly backgated, to control their elec- 
tron density) InAs nanowires. If only SOI is considered, 
the spin degeneracy at the T point (k = 0) is preserved 
because of the time-reversal symmetry. However, this 



degeneracy can be lifted by external magnetic field (fa- 
cilitated in InAs by a large g factor of electrons). When 
the chemical potential is set in the corresponding gap at 
the T point, gapless states which propagate in the oppo- 
site directions with almost opposite spins can be realized 
at the Fermi points. Note that such a system can closely 
resemble the helical edge state of the QSHI [14] . We con- 
sider s-wave superconductor connected to a pair of such 
semiconductor wires in the regime where two CP elec- 
trons split into different wires, in the presence of electron- 
electron repulsion. Effective spin-quantization axes for 
injecting left- and right-moving electrons into the Fermi 
points of the two wires are tilted — in one wire relative 
to the other — by their geometric misalignment. Such tilt 
affects the current cross correlations in the wires in the 
way that is similar to a tunable breaking of the inversion 
symmetry discussed in Rcf. 13J . 

At temperatures and voltage bias between the super- 
conductor and the wires that are smaller than the super- 
conductor gap A, single-particle injection into the wires 
is suppressed. In this regime, transport is dominated by 
the CP tunneling. This process, however, is exponen- 
tially suppressed if the distance between the wires ex- 
ceeds the coherence length of a CP and algebraically on 
the scale of the Fermi wavelength in the superconductor 
(depending sensitively on its spatial dimensionality) [5], 
posing a potentially serious constraint on the interwire 
separation. Very importantly, furthermore, if the applied 
voltage and temperature are smaller than A, the para- 
sitic tunneling of two CP electrons into the same wire 
is suppressed with a power law that is governed by the 
LL correlations 8J. In this work, we thus focus on the 
regime where a CP splits ejecting electrons into the dif- 
ferent wires. There is a time lag of ~ A^ 1 between such 
two tunneling events, the longer it is the weaker the LL 
suppression of the same-wire CP tunneling. However, 
when two electrons are forced to split and enter different 
wires at low energies, the leading-order tunneling rates 
are independent of this time delay (neglecting any inter- 
wire interactions) [8] . Therefore, we consider a simplified 
model with equal-time CP injection of two electrons into 
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FIG. 1. Single-particle electron dispersion with Rashba and 
Dresselhaus SOI. Zeeman splitting 2£ is induced at k = by 
a magnetic field in the z direction, and the chemical potential 
is set in this gap. One-dimensional effective theory is then 
linearized near ±fcp, which define respectively the right- and 
left-moving electron branches. 

two different wires [7J. Note also that electron-electron 
interaction enhances the gap A substantially when the 
chemical potential is tuned appropriately |15j . 

For the wires, Rashba and Dresseslhaus SOI in com- 
bination with the Zeeman splitting are considered. Lat- 
eral confinement in the wire governs subbands, of which 
we suppose (at sufficiently low temperature and appro- 
priate backgate bias) only the lowest to be occupied, 
whose Kramers pairs are split by the lack of both time- 
reversal and inversion symmetries. In this system, the 
one-dimensional effective Hamiltonian for a wire oriented 
along the x axis is given by [TH1 H7J 

H = h 2 k 2 /2m* + aka y + (3ka x - £,& z , (I) 

where m* is the effective mass of electron, a (/3) is the 
strength of the Rashba (Dresselhahus) SOI, and k is the 
electron wave number. The Dresselhaus part is for the 
case when a zinc-blende heterostructure is grown in the 
[001] crystallographic direction, while the wire is oriented 
in the [100] direction [IB]. 2£ = g^sB is the Zeeman 
energy gap at k = 0, with magnetic field B applied along 
the z axis, g is the g factor, and \xb the Bohr magneton. 
<t = (<r x ,a y ,a z ) are Pauli matrices. 

Defining the fc-dependent effective field R(fe) = 
(pk, ak, — £), the Hamiltonian can be written as Hq = 
h 2 k 2 /2m + R(fe) ■ <x, and the eigenspinors are found 
by rotating spinors such that R(fc) • <x|x±(k)) = 
±R(k)\x±(k)), where R{k) = y/k 2 (a 2 + /? 2 )+£ 2 . The 
subscripts +/— here label spin up/down along R. The 
energy eigenstates are thus given by tp±(k) = x±{k)e lkx , 
with energy e±(k) = h 2 k 2 /2m* ± R(k). The upper and 
lower (e + and e_ ) bands are sketched in Fig. [I] When 
the chemical potential /i is set within the gap, we can lin- 
earize the remaining left and right-moving e_ branches 
within a LL picture. This requires (eV, kgT) <C £ and 
electron-electron interactions that are not strong enough 
to hybridize the e± bands. On the other hand, we require 
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FIG. 2. S'-wave superconductor bridging two identical wires. 
The lower wire is rotated by angle 8 with respect the upper 
wire. The superconductor is biased by V with respect to the 
wires. 

the magnetic field to be weak enough on the scale set 
by H c of the superconductor (which can be enhanced in 
mesoscopic structures up to the paramagnetically-limited 
value of order A//z_b [19]). 

Inversion asymmetry between the two wires is intro- 
duced by tilting the lower wire (which is otherwise de- 
fined along the same crystallographic axis), which rotates 
the spin quantization axis at each Fermi point of e_ . The 
upper wire is along the x axis, whereas we suppose the 
lower wire is placed in the xy plane at an angle 9 with 
respect to the x axis, as shown in Fig. [2] (This may 
in practice be realized by growing both wires parallel to 
each other on an unstrained crystal, and then distorting 
the crystal in the xy plane to effectively tilt the wires; 
depending on the interwire separation, a finite 9 may not 
require a large strain, whose additional effect on the SOI 
is neglected.) The SOI in the lower wire is thus given by 

Hso = oik (cos 9<jy — sin^(T x ) + /3k (cos 9a x + sin#<Tj,) . 

This changes the effective fields for the upper (u) and 
lower (d) wires to 

RM (fc) = R(fc, 9 = 0) = (J3k, ak, -f ) , 

RM(fc) = R(M) 
= [k(-asin9 + f3cos9),k(acos9 + j3sm9),-£] . (2) 

The corresponding Fermi-point eigenspinors and 
spin splittings are \x^ u / d \±k F )) = |x^ d ^) and 

R(^ d H±k F ) ee R^. Note that R$ ■ &\ X ^\) = 

— R^l\Xr/i)' an d we wm assume electronic correlations 
are not strong enough to significantly affect these 
Fermi-point spinors. Anticipating tunneling of electrons 
with well-defined spins from the superconductor into the 
Fermi points of our wires, we can effectively decompose 
the fermionic field operators ip^ (<J =t / I) m terms of 
the right/left movers V^/z m the nth wire as |17j 

^ = (xi n) \<j)4 n) + (x^WH { r } ■ (3) 
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The full wire Hamiltonian (JlJ is bosonized [5D] near the 
Fermi points to give an essentially helical (so long as the 
Zeeman term £ is weak) LL [21] : 
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i(M (B) ) a +ff(M B) ) ! 



where ^ B ',9W = (<^>r ± <f>\ )/2 obey commutation re- 
lations = (i7r/2)sgn(a; - y)<W 

parametrize fermionic operators as ip^n oc e ±l ^ r / ! . 

The tunneling Hamiltonian, which describes nonlocal 
injection of the spin singlet CP from an s-wave supercon- 
ductor into the two quantum wires is given by [7] 
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H.c. 



(4) 



In this model, two electrons from a singlet CP split and 
tunnel simultaneously into the upper and lower wires at 
their respective origins. V is the voltage applied be- 
tween the superconductor and the wires, which is set to 
be smaller than A to preclude quasiparticle excitations. 
Expanding spin-dependent operators ipf?9, in terms of the 



chiral modes pertinent to the wires, Eq. (|3|), we can rewire 
tunneling Hamiltonian Q according to 



(d) 



where are the complex-valued expansion coefficients. 

Finally, current-current correlations at the four end 
points of the two wires in Fig. [2] are considered. The 
symmetrized noise spectrum, 



rite*"' ({<*/,(*),<*/,■ (())}>, (5) 



is calculated using Keldysh formalism [T3l [22]. Here 
SIi(t) = Jj(i) — (Ii(t)) are the current fluctuations, i la- 
beling four outgoing channels in the wires (i — 1, up- 
per right; 2 upper left; 3, lower left; and 4, lower right 
branches) . See Fig. [2j These currents are given in the 
bosonic representation by 1^ = {vg/ir)d x 9^ n \ Using 
Eqs. ([3])- ([5]), current correlations are calculated in terms 
of K^, with the final answer (reflecting spin-rotational 
symmetry of a singlet CP) depending only on l-ftT^I 2 , 



\K. 



= (l-RW.R^/2. 



(6) 



Here, r£° 



RjL /i^t , and {K^l 2 can be evaluated 



using Eq. rt2h. R 
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y/k 2 F (a 2 +(3 2 )+e, for n = u, I 



and /i — ±, independent of the orientation of the wire 
or electron chirality. Furthermore, since R 
[ivk 2 F (o? + j3 2 ) cos 9 + £ 2 (identifying 
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we find that |-ftT++| 2 



| 2 and |if+_| 2 



and I = — ), 
\KA 2 . 



Lumping Zeeman and SOI energies into a dimension- 
less parameter A = ^/kp^J ct 2 + /3 2 , we finally arrive at a 
rather simple expression for Eq. ([6]): 

| 2 = (l + A 2 )- 1 (l-^cos6»)/2, 



\K„ 



(7) 



Using these K^, the following expressions for the noise 
spectra are obtained at zero frequency (<j = 0): 
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(8) 



where I = G(eV/e F ) 2 ~>V/(l + \ 2 ), with 7 = (g+g- 1 )/2- 
1, is the total tunneling current from the superconductor 
to the wires (at fc^T <C |eV| ; G oc |r| 2 is proportional 
to the CAR conductance in the absence of LL correla- 
tions) [7JQ2]. This current vanishes in the limit A ^S> 1, 
when both wires become fully spin polarized thus block- 
ing the CP tunneling. Notice that the magnetic field 
did not scramble helical structure of the interwire cross 
correlations, which turn out to be the same [apart from 
the overall suppression by (1 + A 2 )] as in the case of the 
time-reversal symmetric QSHI [13] . This is one of the 
key results of this paper. 

The interwire cross-correlation spectra ([8| are given by 



S±(9, A) cx (l + A 2 )- 1 (l± ff 2 cos6»), 



(9) 



which are modified from those in Ref. [12] only by 
the magnetic- field suppression factor of (1 + A 2 ) -1 . In 
Ref. [13], the angle 9 dependence for the CP injection 
into the helical edge states of a QSHI is due to a tunable 
asymmetry between two edges (induced by a local appli- 
cation of strain or gate voltage to an otherwise inversion- 
symmetric system). Here, 9 dependence comes from the 
mechanical rotation of the lower wire by the angle 9. No- 
tice that the definitions for 5+ and S- are interchanged 
here in comparison to Ref. [13] . This is because the quan- 
tum wires considered here do not have the inversion sym- 
metry of helical edge states on the opposite sides of a 
QSHI strip. Despite this fundamental difference, we can 
clearly see the same structure in the CP noise cross cor- 
relations for both the present quantum-wire system and 
the helical QSHI edges. According to Eq. ([9]), we can ex- 
tract the LL interaction parameter g (which is typically 
g ~ 0.1 — 1 [23] in semiconducting wires) from the inter- 
wire cross correlations: g 2 cos 9 = (S+ — S-)/(S+ + S-). 
While in Ref. [13] the angle 9 is a parameter that may 
not be precisely known, in the present set-up the rota- 
tion angle 9 of the lower wire can be experimentally well 
defined, so that g can be found by measuring S+ and S- 
for an arbitrary value of 9 that is away from ir/2. 

In the discussion so far, we were considering only one 
specific crystallographic orientation of the wires. Namely, 
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the heterostructure growth is in the [001] crystallographic 
direction and each wire is defined (e.g., electrostatically) 
along the [100] direction. However, while the Rashba 
SOI is rotationally invariant around the normal axis, the 
Dresselhaus SOI is sensitive to the wire orientation on 
a crystal's surface JT5J, [53] . Suppose that with the same 
crystal growth direction of [001], the wire is defined at 
an angle 9d from the [100] direction. In this case, the 
Dresselhaus SOI part of the Hamiltonian is given by [TH] 

H D = pk [coB(2B D )a x - sm(26 D )a y ] . 

This crystallographic orientation and the associated 
Hamiltonian are now chosen for the upper wire, with 
our coordinate system still placed (as in Fig. [2]) with the 
x direction collinear with the wire. The corresponding 
effective-field vector is then R^^/c) = [/3k cos(20d), ak — 
/3k sm(20£>), — £]. Since the lower wire is rotated in the 
xy plane by the angle 9 with respect to the upper wire, 
R(') obtained by the corresponding rotation on W- u ' is 
given by RO = [— aksin9 + (3k cos(29d — 6), akcos9 — 
13k sva(29 D -9), -£]. The absolute value of BS U M is mod- 
ified by 9 D : R^ 1 ) = ^Jk 2 \a 2 + /3 2 - 2a/3 sm(29 D )} + £ 2 . 
Both the direction and the magnitude of R("/Q are 
thus modified, affecting in Eq. (pH). We still have 
\K ++ \ 2 = \K—\ 2 and |A + _| 2 = \K-+\ 2 according 
to Eq. (j6j). In fact, the modification of \K^\ 2 can 
be absorbed by redefining A entering Eq. (J7]) as A = 
£/kF\/ce 2 + (3 2 — 2a(3 sin(20o), with all subsequent re- 
lations for the noise spectra unmodified. In particu- 
lar, apart from the modified geometric spin factor A, 
which suppresses the overall strength of the CAR, S± in 
Eq. @ remain the same. This means we can choose any 
wire orientation on the crystal surface without altering 
the essence of the noise cross correlations. One special 
point is 9d = tt/4 when a — (3 (or 9u = —ir/A when 
a = —(3), corresponding to the "persistent spin helix" 
p~6|. 125] , where A blows up and the CAR is fully blocked 
(reflecting exact cancellation of the SOI terms). 

Let us also comment on a possible triplet pairing of 
the injected electrons, e.g., if the two terms in the tun- 
neling Hamiltonian Q acquire a relative phase differ- 
ence: e l5 / 2 ibi u " > tb\ d * > — e^^tbtib^ . We can rewrite it 



is/2)(4 u) 4 d) 



4- 



i(V2)(^Vf 



'). The corresponding modification of |AT M „| 2 
m Eq. can be accounted for by the replacement 
9^9 — 6, with the same S shift of 9 appearing in the 
noise expressions. Interestingly, the phase difference in 
the tunneling terms has the same effect on the current 
correlations as a mechanical rotation of the wires. Such 
triplet component in tunneling can be effectively induced 
by tunneling away from the Fermi points at finite temper- 
ature and/or voltage, and artificially enhanced in more 
complex tunneling setups [25] . 

Another concern to be mentioned is that, if the super- 
conductor is in a slab shape, the perpendicular critical 



field is reduced. This issue can be mitigated by apply- 
ing in-plane magnetic field. For the case of a thin-film 
superconductor, the critical field is further enhanced (up 
to its paramagnetic limit |19j ) when the magnetic pene- 
tration depth is greater than its thickness. R in Eq. |2]) 
needs to be modified accordingly. Since the magnetic- 
field and SOI contributions to R are not perpendicular 
to each other any more, the resulting energy band is not 



symmetric as in Fig. 1 In turn 



in Eq. 



(|7|) and 

the formula in Eq. (|8| acquire some corrections. In the 
limit of A <C 1, the corrections are small, however, and 
we recover the same noise behavior as in Eq. @. In 
the strong-field limit, A > 1, on the other hand, a more 
careful analysis would be warranted. 

Now let us return to Eq. ^ to see the feasability of this 
theory in an experiment. A very large magnetic splitting 
(on the scale of the SOI) in the wires, A 3> 1, blocks An- 
dreev reflection |27j . when the Fermi level is inside the 
T-point gap. The SOI is large in the InAs-based het- 
erostructures and wires, where the Rashba parameter is 
a < 10~ n eVm (being tunable by electrostatic gating) 
[2"8] . (3 <C a, and g factor is ps 15. For electron densities 
in the range of 10 — 100 /im -1 , this gives for the mag- 
netic field B ~ 0.1-1 T corresponding to A ~ 1. Both a 
and g factor can be considerably lower (both up to two 
orders of magnitude) in InGaAs-based heterostructures, 
which can make also a ~ (3 [24] . while the correspond- 
ing magnetic field range remains roughly the same. This 
gives us a favorable operational bound for the magnetic 
field, which opens the T-point gap without compromis- 
ing the strength of the CAR, while also not exceeding the 
paramagnetically- limited critical field (with T c > 1 K). 
Taking everything into account, this means the experi- 
ments can be done at temperatures close to 1 K. 
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